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OO ■ 

' We introduce, in a systematic way, a set of generating functions that solve all the 

different combinatorial problems that crop up in the study of black hole entropy in 
Loop Quantum Gravity. Specifically we give generating functions for: The different 



o 

(N 



sources of degeneracy related to the spectrum of the area operator, the solutions 



to the projection constraint, and the black hole degeneracy spectrum. Our meth- 



< 

O ' ods are capable of handling the different countings proposed and discussed in the 

literature. The generating functions presented here provide the appropriate start- 



ing point to extend the results already obtained for microscopic black holes to the 



Q^, macroscopic regime -in particular those concerning the area law and the appearance 

of an effectively equidistant area spectrum. 
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^r- ■ I. INTRODUCTION 

'NT 

^ ■ The study of the black hole degeneracy spectrum in Loop Quantum Gravity (LQG) has 

OO ! provided important support for the formalism. The confirmation that the expected behavior 

^ ' for the entropy as a function of area is obtained is one of the main physical achievements 
claimed in this framework In addition to the early successes in this respect there has 
^ ■ been an important series of results in the recent past related to this problem. In particular, 
^ ! the studies carried out in [2I, ^ have unearthed a very rich and unexpected behavior of 
the black hole degeneracy spectrum predicted by LQG. For microscopic black holes these 
papers show that, in addition to the exponential growth compatible with the Bekenstein- 
Hawking area law, an effective equi-spacing of the spectrum (with a period approximately 
proportional to log 3) is present. If this feature survives in the macroscopic limit it would 
be a very interesting consequence of LQG because, despite the unevenness of the spectrum 
of the area operator, an emergent, effective, regular spacing would be predicted. This can 
be seen as an additional consistency check for the formalism because such a behavior is 
expected on general grounds 

In order to extend the existing microscopic results to macroscopic areas one has to find 
a way to reach this asymptotic limit without spoiling the content of the theory with un- 
controllable approximations. This is very much in the spirit of Mathematical Statistical 
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Mechanics and Combinatorics. As it is usually done in that framework, the process of 
reaching the asymptotic (thermodynamical) limit requires a number of steps. The first one 
is casting the problem at hand -in this case the counting of the relevant black hole micro- 
states- in such a way that the intimate mathematical nature of the model is captured. This 
has already been done in for the black hole entropy problem in LQG. The success in this 
first step can be judged by trying to carry out the second one: Obtaining suitable generating 
functions for the combinatorial problems involved in the counting of states. This is the 
purpose of this paper. A third -and final- problem, that has to be tackled immediately 
after the one considered here, is to get appropriate asymptotic expansions capturing the 
macroscopic behavior of black hole entropy as predicted by LQG. This may well be the 
hardest step due to its analytic nature. Almost certainly it will require mathematical tools 
different from the number-theoretical and combinatorial methods used to complete the first 
two parts of the program described above. 

The paper is organized as follows. After this introduction we review in section |ll] the 
algorithm proposed in to cornpute the black hole degeneracy. The notation used in the 
paper will closely follow that of |5| . Section [ITIl is devoted to the obtention of the generating 
functions counting the number of solutions to the linear diophantine equations needed to 
describe the degeneracy of the area operator. Section [IV] deals with the generating functions 
giving the full black hole degeneracy spectrum for the different versions of the projection 
constraint that appear in the literature. Finally we end with some conclusions and comments 
in section IVl 

All the results presented here refer to the isolated horizons that are used to model black 
holes in LQG. When we talk about black hole properties in the following we refer, in fact, 
to the isolated horizons representing them. 

II. CHARACTERIZATION OF THE AREA SPECTRUM 

The black hole area in LQG is given by eigenvalues A of the area operator of the form 

N 

/l = 87r74^ VjKjTTi), (2.1) 

7=1 

where 7 is the Immirzi parameter and £p is the Planck length. Here the labels ji G N/2 are 
half-integers associated to the edges of a certain spin network state. They pierce the isolated 
horizon representing the black hole at a finite set of points called punctures [T]. In the 
following we will choose units such that Airjij:, = 1. Horizon quantum states are further 
characterized by an extra label mi that can be interpreted as a spin component. Depending 
on the horizon topology these labels are restricted to satisfy certain constraints that we will 
discuss later. 

The real numbers belonging to the spectrum of the area operator have been character- 
ized in [5|. An obvious, but important, comment is that these numbers must be linear 
combinations of square roots of square-free numbers (SRSFN) pi with non-negative integer 
coefficients g^. In order to check if a number a = Yll=i li^/vl belongs to the area spectrum 
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there must exist jj := ki/2, kj ^ N, such that 

J2 Viki + iy-l = nkVik + iy-l = J2 ^^VPi = «• (2-2) 

7=1 k=l i=l 

Here the denote the number of punctures corresponding to edges carrying spin k/2; hence 
the sum ni + ■ ■ ■ n^^^^ = is just the total number of punctures. Notice that we can always 
write ^^{k^^T}^^^ as the product of an integer times the square root of a square-free 
positive integer number (SRSFN) by using its prime factor decomposition. Equation fl2.2p 
is solved in two steps: First we must identify the allowed spins k/2 such that + 1)^ — 1 
is an integer multiple of some ^/pi, and then determine the value of ra^ that tells us how 
many times each of them appears. In order to deal with the first problem we must solve the 
Pell equations associated to each of the SRSFN's in the r.h.s. of (12. 2p . i.e. 

Vik + 1)^-1 = ^ik + 1)2 - p,y^ = 1, (2.3) 

with ?/ G N. We will label the solutions as {{k^^, y,^) : m G N}, where the index i refers to 
the square-free numbers in each of the Pell equations (see, for instance, for details on the 
Pell equation). Once these numbers are known the can be found by solving the system 
of r-uncoupled jll, linear, diophantine equations 

oo 

Yylnnki^ = 1h i = l,...,r. (2.4) 

m=l 

Notice that, once the are fixed, only a finite number of spins kl^/2, m = 1, . . . , Mj, come 
into play in the equations (12.41) . 

It may happen that some of these equations admit no solutions. In this case Y7i=i 1iy/Pi 
does not belong to the area spectrum. On the other hand, if they do admit solutions, the 
Si=i 1iy/Pi belong to the spectrum of the area operator, the numbers fc^ tell us the spins 
involved, and the n^i^ count the number of times that the edges labeled by the spin kl^/2 
pierce the horizon. 

Let us denote by 5*., z = 1, . . . , r, the set built from the solutions to the i-th diophantine 
equation appearing in (12. 4p as 

m=l 

The elements in these sets are combined in the Cartesian product iS^ = X 1=1 '5*. to give 
all the solutions to the system (12. 4p . The set Sa contains all the spin configurations s = 
{si,S2, . . . ,Sr) G Sa defined by the area a = J2l=i QiVPi- 

Once we have these configurations the black hole degeneracy spectrum is obtained as 

where the sum J2{knk)esi product n(fcnfe)esi extended to the elements {k,nk) of the 
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i-th component of the spin configuration s. The factor P{s) = P{si, . . . , s^) is introduced 
to take into account the projection constraint. The different choices for P will be discussed 
in the following sections. 



III. DIOPHANTINE EQUATIONS 

The purpose of this section is twofold. We will first introduce a generating function 
giving the number of solutions for a collection of uncoupled diophantine equations of the 
form given by (12. 4p . Afterwards we will modify this generating function in order to get the 
reordering degenerations given by the sum of multinomial coefficients in (12.51) obtained by 
taking P{s) = 1 for every configuration. 

The generating functions that we will discuss in the following are written in terms of 
variables that we will denote as Xi with i in one-to-one correspondence with the square- 
free numbers Pi. They involve numerical constants that are obtained from the solutions 
to the Pell equations for each pf, in particular the numbers fc^ and introduced in the 
previous section. For each of the diophantine equations given by (12.41) the generating function 
counting the number of its solutions can be found in any text book on Discrete Mathematics 
or Combinatorics (see, for example, [3]). It has the following simple form 

oo ^ 

^*^°'(^^) = n ; ^- (3-1) 

m=l V-^ -^i J 

The coefficient of xf in the Taylor expansion of (13.11) around Xi = gives the number of non- 
negative solutions to the corresponding diophantine equation (12.41) . Notice that, although 
we are writing an infinite product, in every case we only need a finite number of yl^ (those 
smaller or equal to g^) in order to determine the required coefficient. 

The total number of solutions for a system of such uncoupled diophantine equations is 
just given by the product of the individual generating functions 

oo oo oo , 

G*-\xux„ . . .) = n Gr°\x.) = n n ^^^r- (^.2) 

1=1 i=l m=l l-L -^i J 

Notice again that for a fixed value of the area only a finite number of square-free pi will be 
involved and, hence, the infinite product Hi^i i^, in fact, finite. 

The generating function G*^"^ just computes Xlse^a-'- each allowed value of the area 
(and it gives zero if iS^ = 0). The coefficient of the term xf'--- x^"" in the Taylor expansion of 
Q#so\ jg number of solutions to the system of simultaneous diophantine equations (12. 4p 
and, hence, it coincides with X^se^a-*-' ^^"^ "^^ want to modify (13.21) in such a way that we 
obtain a generating function for the sum of multinomial numbers 



E 



nr T-r I ■ I'J-'Jj 

i=i ll(fc,nfc)esi ""fc- 



A simple way to do it is following a two step approach: First we modify (13.11) to introduce 
the product of factorials in the denominator of (13.31) in front of each term of its Taylor 
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expansion. This can easily be done by considering 

OD OO 



X. 

i=l m=l 



We still have to introduce the factorial appearing in the numerator of (13.31) . This can be 
done by manipulating the previous expression in the following formal way. Let us take 



e"^ exp (Au; Z Z xf^ j dA 

i=l m=l 



G'^"^(^;xi,X2,... 

and consider G^^^{uj\ xi, X2, . . .) for uj = 1. It can be readily seen that 

/ OO OO ^ '^^ 

G\x,, X2, . . . ) = ^^^'^(l; xi, a;2, . . .) = 1 - 5^ 5^ 



i=l m=l 



has the required form. This is a consequence of the following simple formal argument: If 
/(^) = ^'^=o^nX"' then the function g^u), whose Taylor coefficients are n!a„, is given in 
terms of / by 



OO 

-A . 



/OO / OO \ OO 

e-^ J2 anA"cj" dA = ^ nla^cj" . 
\n=0 / n=0 







IV. GENERATING FUNCTION FOR THE BLACK HOLE DEGENERACY 

SPECTRUM 

Let us consider now other choices for P in (12. 5p . Some of them have a direct physical 
meaning whereas others allow us to discuss other possible projection constraints similar in 
form to the standard ones. 



A. Toroidal black holes 

This case corresponds to considering 

r 

pAs) = ll n (k+ir. 

This choice describes a situation in which the third spin components m/ are unconstrained 
and can take any of the kj + 1 possible values independently of each other. In the literature 
this is expressed by saying that no proiection constraint is involved ^ and it can be shown 
that it describes toroidal black holes |8|. The relevant generating function is 

(OO OO \ 

i=l m=l / 
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The coefficient of term xf^ ■ ■ ■ x^"" in the Taylor expansion of the previous expression is the 
total degeneracy (12.51) of a toroidal horizon with area given by qi^/pi + ■ ■ ■ + (ir^fPr- The 
function G^i can be obtained as before in two steps. Consider first 

oo oo 
i=l m=l 

that produces the required inverse factorial terms and also the product ni=i \\.{k rn,)&sS^ + 
I)"'-', and then introduce the factorial term in the numerator of (12.51) by using the same 
formal trick described at the end of section IIIII As in previous instances the formal infinite 
products and sums in (14. ip are, in practice, finite because only a finite number of square-free 
integers are involved for a fixed area value. This means that the generating function can 
be considered, for concrete computations, as a rational function with a finite number of 
variables. 



B. Spherical black holes 

In the case where we have spherical symmetry the so called •projection constraint 



N 



"^mi = (4.2) 



7=1 



must be satisfied by the spin components m/. The accepted view in LQG j^, that we 
will refer to as the DLM counting, is that each m/ is further constrained to satisfy m/ G 
|— fcf /2, ki/2}. There are other proposals in the literature, in particular the GM counting of 
lot, where a different prescription mj G {—ki/2, —ki/2 + 1, . . . , kj/2 — 1, ki/2} is suggested. 



From a purely combinatorial point of view both can be treated in a very similar way so in 
the following we will give generating functions for both approaches. 

The new ingredient that we need in order to take into account the projection constraint is 
a suitable way to count the number of solutions to (14.21) . This can be done in a straightfor- 
ward way. For the standard DLM counting, once the values of kj at the punctures are given, 
the number of solutions to the projection constraint is the constant term in the Laurent 
expansion of 



AT 



7=1 

whereas for the GM counting the number of solutions to the projection constraint is the 
constant term in the Laurent expansion of 

^ ^ ^kj+i ^~ki-i 

7=1 a=0 7=1 



The generating function in these cases can be easily obtained from the toroidal one (14. ip by 
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taking now 



= n n (^'+ 



k r 



n?"(^.^) = n n =11 n 

j=l (fc,nj;)esi a=0 1=1 (k,nk)&Si 



TT TT / Z — Z \ "fc 



fc+1 ^-fc-1 



z — z ^ 



In view of the structure of the P^2^^^{s-, z) and P^^^{s, z) terms we can get the desired gener- 
ating function by substituting the (/c^ + 1) term in (14. ip for (z^'" + z~^"^) or ^ X]a=o 2;'^""^"^ 
respectively. This way we obtain 

G^r{z,x„X2,...) = ' (4-3) 



j=l m=l 

-1 



00 00 



= i-EElE^'^^^'i^'M ■ (4-4) 

\ i=l m=l a=Q J 

These functions have an extra auxihary argument z that is not present in (14.11) . The coeffi- 
cient of the term z'^x^^ ' ' ' tells us the value of the sum (12.51) with a projection constraint 
given by the condition 

m/ = n. 

I 

Notice that, at variance with the cases discussed in the previous sections, the exponents of 
z can be negative and, hence, the expansions that we have to use are Laurent series in z. 
The choice n = corresponds to the spherical black holes. 



V. CONCLUSIONS AND COMMENTS 

We have given a collection of generating functions for a series of combinatorial problems 
related to the description of the black hole degeneracy spectrum in Loop Quantum Gravity. 
The coefficients of their power series expansions give us the exact solution to the counting 
problems that we want to solve. In particular, the generating functions (14. 3 p and (14. 4p give 
us the spherical black hole degeneracy spectrum for the different countings considered here, 
whereas (14. ip gives the one correspondingto the toroidal case. For horizons of higher genus 
it is expected that similar formulas hold [Sj. 

We want to end with some comments. The first is that, despite the apparent infinite 
number of terms involved in the different sums and products appearing in the paper, for a 
given value of area only finite numbers of variables and terms are needed. It is only the fact 
that the diophantine equations that we need to solve have an effective number of variables 
that depends on the area, that forces us to introduce a formally infinite number of them. 

To convince the reader of the power of this generating function techniques we give here 
a concrete numerical example: For an area a = 40-\/2 + 40-\/3 the number of possible 
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configurations can be computed by the considering the generating function 

^ (1 - X\){1 - Xl^){l - X2)(l - 4)(1 - xl^) 

and extracting the coefficient of the term xf^xf^ which has a value of 84. The total degeneracy 
(in the DLM counting) is obtained by taking the generating function 

G°2^{z, xi, X2) = 

1 

~ 1 - (^2 + z-^)xl - (^16 + ;2-16)a;12 _ + ^-1)3,2 - (^6 + ^-6)a;4 _ ^^25 + ^-25)3.15 " 

The value of the black hole degeneracy e'^('^°^+'^°^) is given by the coefficient of the z'^xfxf 
in the power series expansion of G^2^{z, xi, X2). This is 

gS(40v^+40v^) ^ 991809938488860909241077458398212. 

The second comment is that once we have exact closed-form expressions for the black 
hole degeneracies we can ask ourselves about their asymptotic limit and hence extract con- 
clusions for macroscopic black holes. It is very important to realize that without such exact 
and closed-form expressions the problem of extracting all the relevant information in the 
macroscopic limit is very hard and some important features may actually be missed if coarse 
and difficult-to-control approximations are used. An important feature of the black hole de- 
generacy spectrum that one would wish to recover in the macroscopic limit is the effective 
equi-spaced area spectrum found in [2|. In our opinion if such behavior is present it would be 
very strong evidence that LQG provides an accurate description of quantum gravity with the 
right semiclassical limit. We hope that the asymptotic analysis of the generating functions 
given above will uncover this type of behavior. 



Acknowledgments 

We want to thank I. Agullo, J. Di'az-Polo, and E. F. Borja for their comments, encourage- 
ment, and also for providing us with concrete numbers to cross-check the results presented 
here. This work is supported by the Spanish MEG grant FIS2005-05736-C03-02. 



[1] A. Ashtekar, J. Baez, A. Corichi, and K. Krasnov, Phys. Rev. Lett. 80, 904 (1998). A. 

Ashtekar, J. Baez, and K. Krasnov, Adv. Theor. Math. Phys. 4, 1 (2000). 
[2] A. Corichi, J. Diaz-Polo, and E. Fernandez-Borja, Phys. Rev. Lett. 98, 181301 (2007). A. 

Corichi, J. Diaz-Polo, and E. Fernandez-Borja, Class. Quant. Grav. 24, 243 (2007). 
[3] I. Agullo, J. Diaz-Polo, and E. Fernandez-Borja. Black hole state degeneracy in Loop Quantum 

Gravity. arXiv: 0802.3188 
[4] J. D. Bekenstein, Lett. Nuovo Cimento 11, 467 (1974). V. F. Mukhanov, JETP Letters 44, 63 

(1986). J. D. Bekenstein. Classical Properties, Thermodynamics and Heuristic Quantization, 

in Cosmology and Gravitation, M. Novello, ed. (Atlantisciences, France 2000), pp. 1-85. 



9 



[5] I. Agullo, J. F. Barbero G., E. Fernandez-Borja, J. Diaz-Polo, and E. J. S. Villasenor. Black 

hole state counting in Loop Quantum Gravity: A number theoretical approach. Phys. Rev. 

Lett, (to appear), arXiv: 0802.4077. 
[6] D. M. Burton. Elementary Number Theory. McGraw-Hill, New York (2002). 
[7] R. L. Graham, D. E. Knuth, and O. Patashnik. Concrete Mathematics: A Foundation for 

Computer Science. Addison- Wesley Professional; 2nd edition (1994). K. H. Rosen. Discrete 

Mathematics and Its Applications. McGraw-Hill; 6th edition (2006). 
[8] S. Kloster, J. Brannlund, and A. DeBenedictis, Class. Quant. Grav. 25, 065008 (2008). 
[9] M. Domagala and J. Lewandowski, Class. Quant. Grav. 21, 5233 (2004). K. Meissner, Class. 

Quant. Grav. 21, 5245 (2004). 
[10] A. Ghosh and P. Mitra, Phys. Lett. B616, 114 (2005). 



